Three-body Interactions in One Dimension 
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We determine the phase-diagram of a one-dimensional system of hard-core lattice bosons inter- 
acting via repulsive three-body interactions by analytic methods and extensive quantum Monte- 
Carlo simulations. Such three-body interactions can be derived from a microscopic theory for polar 
molecules trapped in an optical lattice. Depending on the strength of the interactions and the par- 
ticle density, we find superfluid and solid phases, the latter appearing at an unconventional filling 
of the lattice and displaying a coexistence of charge-density-wave and bond orders. 



INTRODUCTION 

Quantum many-body systems provide a wealth of fas- 
cinating phenomena in condensed matter physics, includ- 
ing superfluidity in liquid Helium, the fractional quan- 
tum Hall effect, as well as the exotic electronic states in 
the pseudogap regime of cuprate superconductors. While 
these quantum phases emerge from dominant two-body 
interactions, with higher-order many-body interaction 
terms providing only small corrections, an exciting recent 
avenue of research in atomic and molecular physics is to 
engineer systems where higher-order interactions domi- 
nate. In particular, it was recently shown that this goal 
can be achieved for three-body interactions using polar 
molecules [1| . In the present work, we study the most fun- 
damental model Hamiltonian which displays three-body 
interactions in one dimension via quantum Monte-Carlo 
simulations. 

One-dimensional bosonic systems in the strongly corre- 
lated regime have recently been realized with cold atomic 
gases: examples are the superfluid/Mott-insulator quan- 
tum phase transition for atoms trapped in optical lat- 
tices y J 131 1 and the cross-over into the hard-core (Tonks- 
Girardeau) regime (J] . A characteristic feature of hard- 
core bosons in a lattice with additional off-site two-body 
interactions is the appearance of solid phases at half- 
filling n ~ 1/2 with either a charge-density wave (CDW) 
or a bond-ordered (BOW) phase [5[ , (Sl] . 

In contrast, here we study hard-core bosons with 
strong three-body interactions. While the microscopic 
realization of the model with polar molecules gives rise 
to next-nearest-neighbor two-body and three-body inter- 
actions y , the dominant part of the Hamiltonian is 



H 



-J 



E [^I^^+i 



h,b] 



i"i+l 



Ty^nj_injnj+i. (1) 



The first term describes the standard kinetic energy with 
hopping rate J, while the second term, accounts for the 
three-body interaction with strength W; rii = b\bi is the 
density operator with bosonic operators hi and h\ sat- 
isfying the hard-core constraint. A similar model has 



recently been studied in two-dimensions [7|. 




FIG. 1: (Color online) (QMC method: WA) Phase-diagram 
of hard-core lattice bosons with dominant three-body inter- 
actions, Eq. (m, in the grand-canonical ensemble, ^/W vs 
J/W . The solid phase at filling n = 2/3 is characterized by a 
coexistence of CDW and BOW orders. 

We derive the complete quantum phase diagram of 
Hamiltonian ([1]) by means of extensive quantum Monte- 
Carlo simulations. We find the existence of both su- 
perfiuid (SF) and solid phases, see Fig. [TJ However, in 
contrast to systems with two-body interactions, we show 
that the solid phase appearing at the unconventional fill- 
ing n = 2/3 exhibits both CDW and BOW orders. While 
Luttinger liquid theory predicts also instabilities towards 
a solid phase at n = 1/2 and n — 1/3, we show here that 
the system remains superfluid even for strong three-body 
interactions W/ J ^ 1. Solid phases at filling n = 1/2 are 
found by adding weak two-body nearest-neighbor (NN) 
and next-nearest-neighbor (NNN) corrections to Eq. ([T]), 
Vi J2i ^iUi+i and V2 '^^ nini+2 respectively, as naturally 
realized with polar molecules [l|. 



RESULTS 

The ground state phase diagram of the model in Eq. ^ 
is derived in the grandcanonical ensemble by varying the 
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FIG. 2: (Color online) (QMC method: SSE) Lattice filling n = 2/3. (a) Bond-order structure factor S'BOw(27r/3)/L vs W/J. 
Top to bottom: lattice sizes L = 36, 60, 120, 240 and 300; the thermodynamic limit (TDL) is indicated. Dashed line: strong 
coupling perturbative result 4/9( J/W) . (b) Density structure factor Scdw{2tv/3)/L vs W/J. Top to bottom: lattice sizes 
L = 36, 60, 120, 240 and 300; the TDL is indicated. Dashed line: strong coupling result n^/4 = 1/9. (c) Superfiuid density ps 
as a function of W/J, for lattice sizes L = 36, 60 and 120. Inset: Wc/J as a function of 1/ In L, results from the WA. The line 
is a guide to the eye. 



chemical potential fj. at different fixed values of W/J. We 
use two different quantum Monte-Carlo (QMC) meth- 
ods: (i) the stochastic series expansion (SSE) algorithm 
with a generalized directed loop update [8| after a decou- 
pling of the Hamiltonian in trimers for each tree-body 
interaction term, and (ii) a code based on the Worm al- 
gorithm (WA) path integral approach [9|, which allows 
efficient sampling of the many-body path winding num- 
bers in imaginary time and space directions. Although 
the SSE method samples both BOW and CDW orders, 
its efficiency drops for very large values of W/J. The 
WA does not suffer from this issue and is especially use- 
ful in checking the limit W/J ^ 1. Results obtained 
with the two methods are found to be consistent (see be- 
low). The groundstatc properties of the finite systems 
have been probed using temperatures T = 0.6J/L, with 
L the number of lattice sites, which was found sufficiently 
low. 

The phase diagram is determined by two phases, see 
Fig. [T] a superfiuid Luttinger fiquid (LL) phase with al- 
gebraic correlations, surrounding a solid phase at filling 
n ~ 2/3, which appears for dominant three-body interac- 
tions W/J > 3. The solid phase is incompressible, giving 
rise to the characteristic lobe structure in the ii-W phase 
diagram. This incompressible phase is characterized by 
the structure factors Scuw for a charge density wave and 
Sbdw for a bond order at the wave-vector k — 2tt/3 
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with the bond operators Ki — u^^ 

The presence of CDW order can be easily under- 
stood in the limit J = 0, where the ground state at 
filling n = 2/3 is threefold degenerate. In fact, up 
to lattice translations, each groundstate takes the form 
1^) = Ilk^lk^lk+2\^)^ ^nd exhibits CDW order with 
'S'cDw(27r/3) = n L/A. From the unconventional fill- 
ing n = 2/3 of the solid phase, it follows that once 
hopping of particles is allowed the charge density or- 
der implies small hopping correlations, since the position 
of a particle in the solid breaks the symmetry of hop- 
ping to the left or hopping to the right. Using standard 
perturbation theory in the hopping term, we find that 
also a bond order wave appears for finite hopping with 
5Bow(27r/3) == n'^LJ'^/W'^. Monte Carlo results directly 
confirm the coexistence of the two orders. Analogous 
to the case of the Hubbard-model with onsite two-body 
interactions, this incompressible phase at filling n = 2/3 
can be reached by varying the density, which corresponds 
to a mean-field transition, or by keeping the density con- 
stant and by varying the strength of the interactions 
W/J, see Fig. [T] We found that the superfiuid density 
vanishes at the boundaries, indicating that doping of the 
system happens simultaneously to the solid/SF transi- 
tion. In the following we will be mainly interested in 
characterizing the constant-density transition at the tip 
of the lobe. 



S'cDw(fc) = y ^exp[jA:(j - /)] (nj-n/), (2) 
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5Bow(fc) = y^exp[zfc(j-0](/^,A',), (3) 
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Filling n = 2/3 

Our results for the order parameters ^bow and S'cdw at 
fixed density n = 2/3 are shown in Fig.[2]for system sizes 
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FIG. 3: (Color online) (QMC method: WA) Lattice filling n — 1/2. (a) Luttinger Liquid parameter K as a. function of W/J 
with critical value A'c — 0.5 (dashed line) of the superfluid/solid transition. The inset shows the density structure factor 
ScDw(i")/i as a function of the inverse system size 1/L for W/J — 1000. (b) Phase diagram in the presence of an additional 
nearest neighbor interaction Vi in the plane W/J vs Vi/J. The arrow indicates the point of the superfluid/solid transition at 
W/J = 0. 



L = 36,60,120,240 and 300 (top to bottom), together 
with the extrapolated thermodynamic hmit (TDL) be- 
havior. The latter has been obtained based on an ob- 
served linear scaling of 5'bow and Scow in V-^ within 
the solid phase, and is found to be in perfect agreement 
with the strong-coupling results for the order parameters 
valid for W/J ^ 1 (dashed lines). In the thermodynamic 
limit, the two order parameters are found to vanish simul- 
taneously at a critical ratio Wc/J « 2.9, corresponding 
to the solid/liquid transition at the tip of the lobe. 

A more refined determination of the transition point 
can be obtained using bosonization techniques. In the 
weakly interacting regime with W/J <^ 1, the system 
can be mapped onto the sine-Gordon model [l| 
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where $ denotes the charging field and 11 the canoni- 
cal conjugate operator, while K is the LL parameter, 
and 7 = ^/36tt the periodicity of the sine-term. Conse- 
quently, the transition from the LL to the solid phase at 
the tip of the lobe in Fig. [T] appears at the critical value 
Kc = 2/9 and it belongs to the Kosterlitz-Thoulcss (KT) 
universality class. For weak interactions W/J <C 1 the 
behavior of the Luttinger parameter is obtained directly 
from bosonization techniques as K = 1 ~ {2\/3/7r)W/J, 
however this value is strongly renormalized close to the 
liquid/solid quantum phase transition. We compute the 
LL parameter numerically as K = hn^ psn/m by QMC 
methods, with ps and k the superfluid density and com- 
pressibility, respectively. The latter are calculated from 
the statistics of winding numbers {W"^) in imaginary 
time and space directions. In particular, for a square 
system such that Lt ^ L^ ~ i, with Lr = hv /T and 
V = \/ ps/ K.m the sound velocity, ps — mLT{W^)/h'^ [lO|, 



and K = {W^)/LT, respectively, with m = h^ /2J. We 
have performed simulations for L = 60,90, 120, 150 and 
300. In order to precisely locate the critical point we 
employ finite size scaling arguments following from the 
KT renormalization group flow llj. Calling Wc{L) the 
value of W for which K{L) = Kc, the flnite size scaling 
of the transition point is Wc{L) — Wc oc [ln(L)]~^, with 
Wc the transition point in the thermodynamic limit. The 
inset in Fig. [21Ic) shows the finite size scaling of Wc{L), 
which gives the transition point Wc/J = (2.80±0.15), in 
agreement with the discussion above. 



Filling n = 1/2 

In Ref. Ill it is argued that for lattice filling n = 1/2 a 
superfluid/solid transition should occur as a function of 
W/J at a critical LL parameter Kc = 0.5. This insta- 
bility is based on the observation that replacing the den- 
sity operator by its fluctuations around the mean value 
rii = n + drii, the three-body interaction gives rise to a 
nearest neighbor and next-nearest neighbor interaction. 
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In analogy to systems with two-body interactions, one 
would expect a solid phase at half-filling n = 1/2. How- 
ever, we find that the competition between the nearest- 
neighbor interaction of strength 2Wn, which drives an 
instability towards a CDW solid, and the next-nearest 
neighbor interaction of strength Wn, which drives an 
instability towards a BOW solid, removes all instabili- 
ties altogether, see below. The low-temperature phase is 
thus a supcrfiuid. independent of the magnitude of W/J. 
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FIG. 4: (QMC method: SSE) Superfluid fraction ps as a func- 
tion of W/J for fillings n = 1/3, 1/2 and 2/3. Grey and black 
lines correspond to system sizes L = 60 and 120, respectively. 
Due to particle-hole symmetry, the superfluid fractions for 
n = 1/3 and 2/3 are the same for W/J = 0. 



This behavior is a special property of three-body inter- 
actions, and it is in stark contrast to the two-body case, 
where a transition into the sohd phase has been always 
reported [5[. In Fig-EJa) we show our results for K as a 
function of W/ J ^ obtained using the procedure described 
above for system sizes L = 80, 120, 160, 240 and 320. For 
weak three-body interactions W/J < 1, If tends to the 
hard-core value K = 1. We find that for large inter- 
actions W/J ^ 1 the LL parameter saturates at a value 
K = (0.528±0.015) > Kc which is larger than the critical 
value Kc = 0.5, and thus no superfluid/solid transition 
occurs at filling n = 1/2. Consistently, the superfluid 
fraction does not show any appreciable system-size de- 
pendence. Fig. m 

The inset of Fig. OJa) shows the structure factor 
'S'cDw(i') a-s a function of the inverse system size 1/L, for 
the case of W/J = 1000. The structure factor extrapo- 
lates to a value consistent with zero in the termodynamic 
limit, consistent with a superfluid in this limit of large 
three-body interactions. It remains an open question, 
whether there is an analytical result predicting the value 
of the Luttinger parameter for inflnite three-body inter- 
actions in analogy to hard-core bosons with K = 1. 



undergoes a supcrfluid/CDW KT transition at Vi — 2J 
to a two-fold degenerate state with ordering pattern o»o» 
{k = Tr) [13]. 

In Fig. [Sljb) we show how this transition is modifled 
due to the presence of three-body interactions, by pre- 
senting the low-energy phases as a function of W/J and 
Vi/J. We find that the presence of three-body interac- 
tions renormalizes considerably the SF/CDW transition 
by shifting it to lower values of Vi/J < 2. In particular, 
for W/J = 1000 we find V^/ J = (0.75 ± 0.15). 

The presence of NNN interactions V2 > can drive 
the system into a bond-ordered phase. In Ref. |5| it 
is shown that an ensemble of hard-core bosons with 
NNN-interactions (at Vi = 0) enters a k — tt bond- 
ordered phase for V2/J — (2.15 ± 0.10). For even larger 
V2/J > (2.66 ± 0.10), a four-fold degenerate k = tt/2 
CDW phase with ordering pattern o • •o becomes sta- 
bilized. We find that the presence of three-body inter- 
actions significantly renormalizes these transition points. 
For example, at W/J = 10 the SF/BOW transition oc- 
curs for V2/J = (0.9 ± 0.1), and the BOW/CDW transi- 
tion at V2/J = (2.1 ± 0.1). Three-body interactions thus 
widen the extent of the BOW phase in this regime. 



Filling n = 1/3 

We find that at filling n = 1/3 the low-energy phase 
remains a superfluid independently of the strength of 
the three-body interactions. Consistently, the superfluid 
fraction ps in this case is found to be essentially indepen- 
dent of W/J in the range < W/J < 20, as seen in Fig. HI 
We computed the value of the LL parameter for W/J ~ 
1000, and found that K = (0.89 ± 0.01) > K^ = 2/9, 
which confirms that three-body interactions alone do not 
induce a transition in this region. This finding contrasts 
the weak-coupling result Ref. [ij, and exhibits the strong 
breaking of particle-hole symmetry by finite three-body 
interactions, as seen for W y^ from e.g. comparing ps 
for the two fillings n = 1/3 and n = 2/3 in Fig.|4l Adding 
a NN interaction also does not induce any transition at 
n = 1/3. In fact, we find for W/J = 1000 and Vi/J = 20 
that K = (0.48 ± 0.01) > Kc- 



Additional two-body interactions 



CONCLUSIONS 



We finally consider the effects of additional short- 
ranged two-body repulsions, by adding to the model in 
Eq. (1) a term Vi J2i ^i'ni+i + ^2 X^i "i''^i-i-2 with nearest 
neighbor (NN) repulsion Vi, and next-nearest-neighbor 
(NNN) repulsion V2. Such terms are important as they 
constitute the leading correction to the short-range three- 
body repulsions considered thus far, and can induce new 
instabilities of the LL. In fact, it is well known that a half- 
filled system of hard-core bosons with NN repulsions Vi 



We have determined the phase diagram of the fun- 
damental model Hamiltonian for hard-core bosons with 
dominant short-range three-body interactions in one di- 
mension. The latter can be realized with polar molecules 
cooled to the electronic and vibrational groundstate, 
which is well in the reach of current experiments [l3| . The 
one-dimensional nature of the problem opens fascinat- 
ing prospects for studying the dynamics of systems with 
many-particle interactions by using powerful numerical 



techniques such as e.g. tDMRG IJ]. Extensions to two 
dimensions [7| in various lattice geometries hold promises 
in the search for exotic phases, such as e.g. topological 
phases and spin liquids. 
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